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Tj- ; ABSTRACT 

oo 

' Context. Cygnus X-1 (Cyg X-1) is a high mass X-ray binary system, known to be a black hole candidate and one of 

the brightest sources in the X-ray sky, which shows both variability on all timescales and frequent flares. The source 
spends most of the time in a hard spectral state, dominated by a power-law emission, with occasional transitions to 
the soft and intermediate states, where a strong blackbody component emerges. 

Aims. We present the observation of Cyg X-1 in a hard spectral state performed during the AGILE science verification 
phase and observing cycle 1 in hard X-rays (with SuperAGILE) and gamma rays (with the gamma ray imaging detector) 
and lasting for about 160 days with a live time of 6 Ms. 

Methods. We investigated the variability of Cyg X-1 in hard X-rays on different timescales, from ~ 300 s up to one day, 
and we applied different tools of timing analysis, such as the autocorrelation function, the first-order structure function, 
?— ( ' and the Lomb-Scargle periodogram, to our data (from SuperAGILE) and to the simultaneous data in soft X-rays (from 

5^ , RXTE/ASM). We concluded our investigation with a search for emission in the energy range above 100 MeV with the 

maximum likelihood technique. 

Results. In the hard X-ray band, the fiux of Cyg X-1 shows its typical erratic fluctuations on all timescales with variations 
of about a factor of two that do not significantly affect the shape of the energy spectrum. From the first-order structure 
function, we find that the X-ray emission of Cyg X-1 is characterized by antipersistence (anticorrelation in the time 
series, with an increase in the emission likely followed by a decrease), indicative of a negative feedback mechanism at 
work. In the gamma ray data a statistically significant point-like source at the position of Cyg X-1 is not found, and the 
upper limit on the flux is 5 x 10~* ph cm"'^ s~^ over the whole observation (160 days). Finally we compared our upper 
limit in gamma rays with the expectation of various models of the Cyg X-1 emission, both of hadronic and leptonic 
origin, in the GeV - TeV band. 

Conclusions. The time history of Cyg X-1 in the hard X-ray band over 13 months (not continuous) is shown. Different 
analysis tools do not provide fully converging results of the characteristic timescales in the system, suggesting that the 
timescales found in the structure function are not intrinsic to the physics of the source. While Cyg X-1 is not detected 
in gamma rays, our upper limit is a factor of two lower than the EGRET one and is compatible with the extrapolation 
of the fiux measured by COMPTEL in the same spectral state. 
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1. Introduction 

Cygnus X-1 (Cyg X-1) is one of the brightest X-ray sources 
in the sky. It is a binary system eomposed of a compact 
object and the 09.7 lab supergiant star HDE 226868, fill- 
ing 97 % of its Roche Lobe, with a m ass ranging be- 
tween 15 and ^ 30 (see for e xample iGierlinski et al.l 
Il999t ICaballero-Nieves et al.l I2009D . The measurement of 
the mass of the com pact object, with a range between 4.8 
Mq a .nd 14.7 Mp, bvlHerrero et al.l (I1995D and 8.7±0.8 Mq 
from IShaposhnikov fc TitarchukI ( 20071 ) . suggests identifi- 
cation with a black hole. The di stance to th e binar y system 
is measured as 2.1 ± 0.1 kpc bv lZiolkowskil (120051 ). 

A characteristic feature of the black hole binaries (as 
Cyg X-1) in the X- ray band, discussed for example by 
iFrontera et al.l ()2001h between 0.5 and 200 keV, is the exis- 
tence of two well distinct emission states: "low/hard" and 
"high/soft". The typical energy spectrum in the low/hard 
state, in which the source spends most of its time, is 
described well by a power-law {E~^) with photon index 
r ^1.7 and a high-energy cutoff at ^ 150 keV. Instead, 
in the high/soft state the source is characterized by a 
strong blackbody component with kT ~ 0.5 keV and a 
soft power-law tail with F usually ranging between 2 and 3. 
An "intermediate" spectral state also exists, discovered by 
iBelloni et al.l (jl996l) in observations with the Proportional 
Counter Array (PCA) aboard the Rossi X-Ray Timing 
Explorer (RXTE), in which the flux is higher of about a 
factor of two with respect to the low/hard state and the 
spectrum is softer (with a photon index of ~ 2.1 and a 
blackbody component of ~ 0.36 keV temperature). 

It is useful to note that the definition of low/hard and 
high/soft states derives from the observations at soft X- 
rays. When observing in hard X-rays, the condition re- 
verses: the source is a factor of ^2 brighter in the low/hard 
state than in the high/soft one. The typical fiux in the 
low/hard state is 8 x 10^^ erg cm~-^ s^^ in 20 - 40 keV 
and ~ 3 X 10~® erg cm~^ s~^ in 1 - 10 keV. Although these 
definitions may be misleading when applied to the obser- 
vations in hard X-rays reported in this paper, throughout 
the text we use the classification low/hard versus high/soft 
to comply with the classical literature on this source. 

Cyg X-1 is a highly variable source in X-rays. Its vari- 
ability is observ ed on any timescale, from months to sec- 
onds (see e. g. Brocksopp et al.l 119991 : iPottschmidt et all 
[20031: iLing et al.l 11991 " In particular, in the hard X-ray 
range the experiments of the Interplanetary Network de- 
tected seven episodes of giant flaring in the 15 - 300 keV 
energy range, with a duration of 0.9 to 28 ks, peak fiux of 
order of 10~^ erg cm~^ and fiuence ranging from 5 x 
10~^ erg cm~^ to 8 x 10~* erg cm~^ (see iGolenetskii et al.l 
|2003( ). These outbursts were detected during both low/hard 
and high/soft spectral states, and, in general, the giant 
bursting events seem to maintain the spectral parameters 
(and likely the emission mechanisms) of the underlying 
state. 

Re cently, Cyg X-1 has been observed by lAlbert et al.l 
(|2007l) above 150 GeV energy with the Major Atmospheric 
Gamma Imaging Cherenkov (MAGIC) telescope. A TeV 
excess of 4.1 a compatible with a point- like source and spa- 
tially consistent with Cyg X-1 was observed simultaneously 
with a hard X-ray fiare taking place in the low/hard state 



(-1.5 Crab - 1.2 x 10"^ erg cm"^ s"^ - in 20 - 40 keV 
with INTEGRAL ("Malzac et al.l l2008l) and ~1.8 Crab in 
15 - 50 keV of Swift/BAT and ~0.6 Crab in 2 - 12 keV of 
RXTE/ ASM). The TeV excess was detected at the rising 
edge of the hard X-ray peak, one day before its maximum, 
while no variation is found in the soft X-ray emission. The 
source does not show any steady emission in the TeV band 
and the upper limits above — 150 GeV energy at the 95 % 
confidence level reach 2 % of the Crab Nebula fiux. 

Although Cyg X-1 is a well-known source in X-rays, very 
little information is known about its emission in gamma 
rays. During the first three cycles of observation (1991 - 
1994) of the Compton Gamma Ray Observatory (CGRO), 
the source, in hard state, w as detected by COMPT EL only 
between 2 and 5 MeV (see iMcConnefl erall[2000l for de- 
tails). EGRET did not detect the source during that ob- 
servation and the upper limit to the fiux is of the order of 
10 X 10"^ ph cm~^ s~^, posing no need for a high-energy 
cut-off. In the soft state, the spectrum of Cyg X-1 can be 
modelled with a power-law that extends with the same pho- 
ton index (~ 2.5 - 3) beyond 1 MeV and up to about 1 
MeV as detected by COMPTEL (Mc Connell et al.ll2002D . 
Unfortunately in this case, no EGRET measurement is 
available. 

In steady conditions, the radio e mission of Cyg X-1 i s 
stable during low/hard states (s ee iGleissner et alll2004D . 



except for rarely observed flares ([Fender et al.l 20061) . and 



appears to be q uenched below a detectable level during the 
high/soft state (j Brocksopp et al. 1999). A non-thermal ra- 
dio jet, extending up to — 15 x 10^'^ arcsec with an opening 
angle less than 2°, w as detected in VLBA observations by 
IStirling et all (l200lD . 

The X-ray flares and a radio-emitting jet allow the 
cla ssiflcation of Cyg X-1 as a microquasar, as described 
by iMirabel fc RodrlguezI (jl999l) . The high-energy parti- 
cles in the radio-emitting jet ar e likely to produce gamma 
rays as well (see Dubud 20071) . Following the model by 
Zdziarski & Gierlihsk i ni2004D . another possible source of 
gamma rays, especially in the low/hard spectral state, is the 
high-energy tail of the electron distribution in the corona. 
Instead in the high/soft state, the energy spectrum does 
not show an energy cutoff, a s conflrmed by the COMPTEL 
detection up to ~ 10 MeV (jMcConnefl et al.l [200211 and a 
gamma ray emission above thi s energy is expected. 

AGILE ([Tavani et al.ll2008D is the first satellite mission 
sensitive to gamma rays (in 30 MeV - 30 GeV) flown af- 
ter EGRET. It observed Cyg X-1 for - 160 days in four 
different pointings during its first obse rving cycle. Th anks 
to the X-ray monitor Super AGILE (Feroci et al.ll2007l) . the 
source is observed at the same time in the hard X-ray band 
(between 20 and 50 keV) . AGILE can continuously monitor 
the source, with a duty cycle of about 50 %. Other scanning 
instruments, such as RXTE/ ASM or Swift/BAT, observe 
each source many times a day but for a shorter duration, 
and the duty cycle is usually shorter than 10 %. 

In this paper we report the probably longest uninter- 
rupted observation of Cyg X-1 in the hard X-ray and 
gamma ray bands from the AGILE data, complemented 
and extended at lower energy (2 - 12 keV) with the in- 
formation from the public web archiv43[3 of the All Sky 
Monitor (ASM) aboard RXTE. After the summary of the 
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observations (given in Sect. 2) and the description of the 
methods of analysis (Sect. 3), we report the results of the 
data analysis in X-rays (Sect. 4) and gamma rays (Sect. 
5). Finally, in Sect. 6 we discuss our findings and draw our 
conclusions in Sect. 7. 



2. Observations 

AGILE observed the Cygnus region four times during the 
first cycle of observations: from early November until mid 
December 2007 (44 days, hereafter observing block 1 or 
OBI), in the second and third decades of April 2008 (20 
days, hereafter observing block 2 or 0B2), in the second 
half of May and whole June 2008 (51 days, hereafter ob- 
serving block 3 or 0B3), and from mid October until the 
beginning of December 2008 (45 days, hereafter observing 
block 4 or 0B4). The first observing block began during 
the satellite science verification phase, which lasted from 
July until November 2007. 

Each observation in the AGILE pointing plan usually 
lasts from 2 to 4 weeks to favour the photon collection in 
the gamma ray band. The main constraint in the AGILE 
pointing strategy is the angle between the Sun and the so- 
lar panels, which has to be fixed at 90° ± 1°. Because of this 
solar constraint, the pointing of AGILE, hence the position 
of the sources in the field of view, drifts by ~ 1° per day. 
For this reason the off-axis angle (the distance between the 
source position and the satellite boresight) is not constant 
during the observation. The value of the off-axis angle of 
Cyg X-1 during the four observing blocks ranges from 2.4° 
to 32.7°. The exposure of each observing block, as reported 
in the table, is computed as the total time in which the 
source is observed, excluding the duration of the satellite 
passages through the South Atlantic Anomaly (SAA) and 
the source occultation by the Earth. Following these crite- 
ria, the total on-source exposure amounts to 6.1 x 10^ s. 
We show in table [T] the details (start and end date, off-axis 
angle and exposure) of the observing periods in the AGILE 
observation. 



3. Data reduction and analysis 

3.1. X-rays 



SuperAGILE (see lFeroci et al.|[2007l for a description) is a 
coded aperture instrument, and its data analysis software 
is highly customized, as frequently happens for similar ex- 
periments. The data of SuperAGILE are processed using 
the Enhanced Multi sky Imaging (EMI) software package, 
suitable for decoding the ir aages of bright sou rces in mod- 
erately crowded fields (see iFeroci et al.l l2010l for details) . 
In the observation reported in this paper, Cyg X-1 is al- 
ways the brightest source in the field of view, making EMI 
a reliable method of data processing. 

SuperAGILE is a 1-D imager with a field of view 
composed of two orthogonal and overlapping regions, of 
107° X 68° each at zero response. In the central area, of 
68° X 68°, both coordinates are encoded, thus giving twice 
a 1-D imaging, and the sensitivity on axis is of the order 
of 15 mCrab at 5ct in one day with an energy resolution of 
8 keV FWHM. During all the observations, Cyg X-1 is 
in the twice 1-D region, with an off-axis angle ranging from 
2.4° to 32.7° (see table [I]). As stu died in the SuperAGILE 
ground and in-flight calibrations (jEvangelista et ahl 120061 : 



iFeroci et al.ll2008h . the variation of the instrument Point 
Spread Function in this range of off-axis angles does not 
affect significantly the imaging response. 

The output of the EMI software, for all the sources de- 
tected in an image above a significance level of 5cr, is the 
count rate normalized with the instrument effective area 
convolved with a Crab-like spectrum, making it indepen- 
dent of the source position in the field of view. Throughout 
this paper we refer to the normalized counting rate defined 
above, and as a reference, the value for the Crab Nebula 
is 0.15 cts cm~^ s~^. The normalized rates from EMI have 
been calibrated by means of a raster scan of pointings to- 
ward the Crab Nebula during the AGILE science verifica- 
tion phase and also verified with other bright sources using 
the Swift/BAT transient monitor resultf0. The uncertainty 
on the normalized rate values quoted throughout this pa- 
per includes both the s tatistical and systematic components 
fsce Fcr oci et al.l l20lJ. for details). Since the systematic un- 
certainty applies to the absolute value of the flux but not 
to relative values of subsequent measurements in the same 
pointing, in the plots we show only data points with statis- 
tical errors, separately showing the systematic uncertainty 
on the flux calibration. 

We accumulated the "correlated" energy spectrum of 
Cyg X-1, which is the flux of the source extracted from im- 
ages integrated in different energy intervals. The spectrum 
of Cyg X-1 was extracted in bins of 2 keV amplitude, from 
20 to 50 keV, by using the EMI software. Since the source 
is bright, the significance of the source in the images is 
greater than 5a in all the energy intervals. The spectrum is 
accumulated separately in the two 1-D directions and then 
fitted s imultaneously using the software package XSPEC v. 
12.3.0 (jDorman et al.ll2003D . 



3.2. Gamma rays 

We completed the study of Cyg X-1 with the observation by 
the AGILE gamma ray imaging detector (GRID), sensitive 
in the energy band between 30 MeV and 30 GeV. The data 
are anal ysed with the AGILE Standard Analysis Pipeline 
(see e. g. IVercellone et aI1l2008l for a detailed description). 
We accumulated the count map, calculated the exposure 
map, and estimated the map of the gamma rays from the 
Galactic Background only above 100 MeV and for events 
flagged as confirmed gamma rays (filter code=5) after ex- 
cluding the data collected inside the SAA (phasecode=18). 
The region of interest in the study of Cyg X-1 reported 
in this paper has a 40° radius, inside a map of 80° ra- 
dius with 0.5° bin size. We used a particular application of 
the AGILE maximum likelihood (ALIKE) procedure, called 
ALIKEsingle, to search for the detection of sources at a 
specified position. Every run of the ALIKEsingle software 
provides the user with the significance of the maximum like- 
lihood technique applied at the input position and either 
the flux or the 2a upper limit in case a point-like source 
is detected or not, respectively. To sum all the data of the 
whole gamma ray observation, all the maps are centred on 
the same galactic coordinates {I = 88.0°, b ~ —12.0°). 
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4. Results in the X-ray band 

4.1. The hard X-ray variability 

Thanks to the continuous monitoring of SuperAGILE, a 
measure of the variabiUty of Cyg X-1 with the resolution 
of about 3 ks from the standard EMI pipeline is obtained. 
We adopted this resolution, corresponding to an AGILE 
orbit, to guarantee the statistical significance of the source 
detection in each time bin and to favour the joint analysis 
with RXTE/ASM, although the SuperAGILE data allow 
for finer time resolution. A synthetic view of the normal- 
ized counting rate measured by SuperAGILE is given in 
Fig. [3 together with the measure of RXTE/ASM in the 
same time interval. The time bin of ^ 3 ks is used for both 
instruments. As a reference, the counting rate of the Crab 
Nebula is 0.15 cts cm~^ s~^ in SuperAGILE and 75 cts s^^ 
in RXTE/ASM. 
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Fig. 1. Superposition of the SuperAGILE (upper panel) 
and RXTE/ASM (lower panel) fluxes of Cyg X-1 during 
the observation. The error bars are omitted for clarity and 
the RXTE/ASM data are rebinned with 3 ks time bin. The 
vertical dashed lines limit the observing blocks. The nor- 
malized rate of the Crab Nebula is 0.15 cts cm^^ in 
SuperAGILE and 75 cts s'^ in RXTE/ASM. 



integration between MJD 50087 and 54801) is represented 
by grey dots and we superimpose the data simultaneous to 
the AGILE observation (black dots), the points (squares) 
during the soft state between MJD 50230 and 50307 (see 
IZdziarski et al.|[2"002|) and during the interme diate state of 
MJD 52797 - 52801 (crosses), reported by iMalzac et all 
(|2006f ). From the distribution of the black points in the 
colour-colour diagram, it is possible to see that Cyg X-1 
remained in the low/hard spectral state for the whole du- 
ration of the AGILE observation. 
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So-'t Rotio (5 - 12 keV / 1.5 - 3 keV) 

Fig. 2. Colour-colour diagram of Cyg X-1 from the public 
data of RXTE/ASM with the hard ratio (estimated from 
the counting rate in 5 - 12 keV divided by the one in 1.5 
" 3 keV) as a function of the soft ratio (from the fraction 
between th e rate in 3 - 5 keV and the one in 1.5 - 3 keV), 
similarly to iReig et al.l (|2002[ ). The grey dots are the daily 
average of the whole dataset (from MJD 50087 until 54801), 
the black dots are the daily average of the data simultane- 
ous to the AGILE observation, the black squares are the 
data in the soft sta te from MJD 50230 until 50307 (see 
IZdziarski et al.|[2002[ ) and finally the black crosses are the 
dwell averages i n the intermediate state from MJD 52797 
until 52801 (see lMalzac et alll2006l) since this state lasted 
for only five days. 



The continuous and uninterrupted source covering of 
SuperAGILE allows us to detect possible fiaring episodes 
and also to study the fiux with higher time resolution, if 
compared for example to Swift/BAT with its higher sensi- 
tivity in a similar energy band (15 - 50 keV) but a com- 
pletely different pointing strategy with shorter "snapshots" 
at higher sensitivity but lower duty cycle. 

By comparing the rate of Cyg X-1 measured by 
RXTE/ASM during our observation with typical values 
in the hard spectral state (20 - 30 cts s~^) and in the 
transitions to the soft spectral state (more than 80 - 100 
cts s~ ^ on June 1996 (iFrontera et al.lT200l[ ) or September 
2001 (|Pottschmidt et al.l 120031 )1" we derive that the source 
did not reach the soft state during our observation. We 
study with more completeness the spectral state during the 
AGILE observation by means of a colour-colour diagram 
from the pubhc data of RXTE/ASM, plotted in Fig. ^ In 
the figure the whole dataset of RXTE/ASM (with daily 



We select two different timescales to estimate the source 
counting rate in the SuperAGILE band: one satellite orbit, 
with net exposure of ~ 3 ks, and one day, with exposure 
of ~ 40 ks. The normalized rate of Cyg X-1 on orbital 
timescale is plotted in Fig. [31 All the plots have the same 
vertical scale to compare the source variability. The rate 
is converted into units of erg cm^^ using a conversion 
factor derived from the fit of the average energy spectrum. 

In the plots the characteristic slow variations of the 
emission of Cyg X-1 are clearly detected, with a typical 
timescale of 30 - 50 days in which the source flux changes 
up to a factor of ^ 3. Superimposed to this slow modula- 
tion, a faster and erratic variability, with timescale of about 
one day and a maximum flux variation of the order of two, 
can be seen. In the same figures it is also possible to appre- 
ciate the short time flickering of Cyg X-1, on the scale of 
a few hours. Although some periods of rapid fiux increase 
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were detected, for example around MJD 54420, 54605 or 
54648, large flares were not detected during the AGILE ob- 
servation. 

SuperAGILE always stores and transmits all the de- 
tected counts (photon- by- photon operative mode). We 
could then measure the normalized rate of Cyg X-1 ac- 
cumulating images with a time bin of ~ 300 s exposure, a 
trade-off between time and statistical quality of the images. 
We excluded the time intervals corresponding to the SAA 
and to Earth occultation. Usually we cut the SuperAGILE 
data with an energy threshold fixed at 20 keV to avoid the 
systematics due to the variation of the instrument tempera- 
ture. Instead, in the integration time of ~ 300 s the temper- 
ature, and consequently threshold, variation is small, thus 
we can lower the energy threshold down to 17 keV. We show 
for example in Fig.|3]the lightcurve of Cyg X-1 with a reso- 
lution of - 300 s between MJD 54604 and 54605, during a 
flaring episode. The higher time resolution allows the short- 
time variability superimposed to the main flux variation to 
be appreciated. 
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Fig. 4. SuperAGILE lightcurve of the source flux increase 
between MJD 54604 and 54605 in the energy interval 17 - 
50 keV with a bin size of ~ 300 s. The gaps in the lightcurve 
are due to the passage of the satellite through the SAA and 
to the source occultation by the Earth. Only the statistical 
uncertainty is shown on the points in the plots. The average 
systematic uncertainty on the flux calibration is plotted in 
the top right angle. 



The availability of the photon-by-photon data also per- 
mits building of the energy spectra. To check for spectral 
variability possibly associated with the observed flux vari- 
ations, we studied the energy spectrum of Cyg X-1 in three 
flux states: average (from MJD 54610 until 54612), faint 
(between MJD 54597 and 54599), and bright (MJD from 
54761 to 54762). A simple power-law, th e model pegpwrlw 
of XSPEC v.12.3.0 (|Dorman et al.ll2003D . provides a good 
fit to the spectrum in the 20 - 50 keV energy band (see Fig. 
[5]). The results of the spectral fit are summarized in Table 
[H The spectral shape of the source does not show signifi- 
cant variations despite a flux variation of nearly a factor of 
two, except for a marginal indication of a small softening 
at lower flux. 




Fig. 5. Spectrum of Cyg X-1 from the SuperAGILE im- 
ages in the X direction showing the bright (MJD 54761 
- 54762, 75 ks exposure, 0.17 cts cm~^ s~^, red in colour 
version) and faint (MJD 54597 - 54599, 92 ks exposure, 
0.08 cts cm~^ s^^, black in colour version) flux states of 
the source. The parameters of the fit for the images in both 
directions are reported in Table [2J where the uncertainties 
are at a 90 % confidence level and include a 5 % systematic 
error. 



4.2. Correlation between soft and hard X-rays 

We studied the source variability in the SuperAGILE data 
with the sampling time of one day, by accumulating the im- 
ages separately in two energy bands, low energy (LE) from 
20 to 25 keV and high-energy (HE) from 25 to 50 keV, with 
a net exposure of ~ 40 ks for each image. A hardness ratio 
is computed as the rate in HE divided by the one in LE. 
The plots of the normalized rate and hardness ratio are 
shown in Fig. [51 We find that, during the observing blocks 
1, 2 and 4, the hardness ratio is roughly constant despite 
rate variations by roughly a factor of ^ 2, with a few ex- 
ceptions; for example, between MJD 54441 and 54443, the 
hardness ratio increases in coincidence with an enhance- 
ment of the counting rate after a local minimum, owing to 
a flux reduction of about 50 % on MJD 54442. In the ob- 
serving period 3, the hardness ratio shows a general linear 
increase, with short-term fluctuations superimposed, cor- 
responding to a similar trend in the rate. The comparison 
with the data in 2 - 12 keV (see Fig.[T]) shows that the soft 
X-ray flux decreases at the same time thus confirming the 
higher hardness. 

We investigated the contribution of the emission mecha- 
nism at different energy ranges by studying the correlation 
between the soft (from the RXTE/ASM public archive) and 
hard (SuperAGILE observation) X-ray bands. We removed 
the dips simultaneous to th e inferior coniunction of the bi- 
nary system (observed by 'Maso n et al.l 119741 : iHolt et aD 
1979; Pricdhorsky et al. 1995), with the stellar companion 
passing between the black hole and the observer along the 
line of sig ht, using the orbita l solutions for HDE 226868 re- 
ported bv lLasala et al.l()1998D , with a period of 5.5998 days. 
We then rebinned the lightcurves of both instruments on a 
common timescale with the same bin size of six hours, in 
order to increase the statistic of the measurement and im- 
prove the chance of having at least one ASM measurememt 
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in the bin. From the superposition of the lightcurves and 
the hardness ratio for each observing block, shown in Fig. 
[71 we can see that there are intervals of time of increased 
fluctuation in the hardness ratio. For example in the OBI, 
the hardness ratio is almost constant between MJD 54406 
and 54427, and its variability increases until MJD 54442. 
Simultaneous to the increased variability, we can see that 
the hard X-ray flux shows a flaring episode of about five 
days duration reaching ~ 800 mCrab around MJD 54438, 
while the soft X-ray flux has a decreasing trend. A simi- 
lar behaviour is shown in the 0B4, with an increase in the 
hardness ratio between MJD 54785 and 54793, correspond- 
ing to a bump in hard X-rays, peaking around 500 mCrab 
on MJD 54788, while the soft X-ray emission is almost con- 
stant at ~ 250 mCrab level. At the end of the 0B3, the flux 
measured by Super AGILE increases up to ^ 1500 mCrab 
while the measure by RXTE/ASM remains between 200 
and 400 mCrab, but in this period the source is at an off- 
axis angle of about 32°, where the flux estimation is less 
accurate. Unfortunately, the RXTE/ASM data are sparse 
in these time intervals and consequently we cannot draw 
any serious conclusion. To verify if we can find any varia- 
tion in the source photon index, we extracted the energy 
spectra of these episodes, but we did not find any signifi- 
cant differences from the ordinary values. For example, the 
photon index of the power-law that fits the spectrum ac- 
cumulated in the 0B4 between MJD 54787 and 54789 is 
consistent with the values reported in Table [2j 

To investigate possible correlations between the source 
emission in soft and hard X-rays, we show in Fig.|8]the scat- 
ter plot of the fluxes of SuperAGILE and RXTE/ASM, each 
one converted into Crab units for clarity. A certain degree 
of correlation appears in the plot, although the value of the 
correlation coefficient is 0.03 and does not allow establish- 
ment of a formally significant correlation, even with more 
than 600 degrees of freedom. In the lower lefthand region 
of the plot (between ~ 150 and ^ 250 mCrab values for 
the ASM rate), we can see the data corresponding to the 
intervals in which the source emission is increased in the 
hard X-ray range and is almost constant in soft X-rays. We 
also computed the cross-correlation of the RXTE/ASM and 
SuperAGILE lightcurves (both accumulated with six hours 
resolution), and we did not find any time lag greater than 
the bin size of six hours. The degree of correlation at zero 
time lag is higher in the OBI (0.73) and 0B2 (0.69) than 
in the 0B3 (0.54) and 0B4 (0.55). 

4.3. Analysis of the structure function 

The SuperAGILE data are not well-suited to timing analy- 
sis of Cyg X-1 using the standard tools, e.g. the fast fourier 
transform. In fact, the measured rate is a few counts per 
second, depending on the position of the source in the field 
of view, and the occultation of the source by the Earth pro- 
duces a discontinuity (and consequently a windowing effect) 
on the time series. 

As explained for example bv lHughes et al.l (|1992[ ). a use- 
ful tool for studying the time variability of a physical quan- 
tity with such sampling characteristics is the first-order 
structure function, simply referred to as "structure func- 
tion" or SF. This function is commonly used in the anal- 
ysis of time series (see e.g. iRutm an 197 ^, has been intro- 
duced in astronomy bv ISimonetti et alT l 1985D mainly for 
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Fig. 8. Scatter plot of the SuperAGILE (20 - 50 keV) vs 
RXTE/ASM (2 - 12 keV) flux, both separately converted 
into Crab units and rebinned on a common timescale with 
a bin size of six hours. The typical uncertainty is plotted 
for clarity only in the lower left corner. 

temporal gaps in the data. The insensitivity of the struc- 
ture function to gaps in the data is extremely important for 
satellites in low earth orbit, such as AGILE, whose obser- 
vations are interrupted by the occultation of the source by 
the Earth and by the passage of the satellite through the 
South Atlantic Anomaly. 

The first-order SF is defined by ISimonetti eraD (|1985[) 
as Dj{T) = {[f{t + t) — fit)]"^) and in our case it can be 
estimated as 



(1) 



the study of active galactic nuclei, and is also insensitive to 



where the brackets () indicate the average, t is the time 
lag, f{t + r) and f{t) represent the fiux at time t + t and 
i, respectively, and N{t) is the number of pairs considered 
in the average. 

A plateau at small time lags (often referred to as "hor- 
izontal branch" ) is common of the SFs of the active galac- 
tic nuclei in optical (see for example iPaltani et al.l 1997f) 
and gamma rays (see iNandikotkur et al.lll997r and is gen- 
erally produced by the white noise (uncorrelated fiux) of 
the source emission, yielding a constant value independent 
of the time lag. A monotonic region is found at increas- 
ing lags, where the SF approxi mately follows a power-law , 
Djir) oc t'^" (see for example INandikotkur et al.l[l997l ). 
The parameter H is also known as the Hurst exponent and 
is used in describing non-stationary processes. Systems with 
0.5 < H < 1 are persistent: an increasing trend in the 
past is probably followed by a similar increasing trend in 
the future. Conversely, if < H < 0.5 the process shows 
antipersistence and an increase in the past is probably fol- 
lowed by a decrease in the future. Finally, the minima in the 
structure function indicat e typical timescales of the system 
(INandikotkur et al.l[T997l) . 

We computed the structure function of the normalized 
rate of Cyg X-1, extracted separately in the four observ- 
ing blocks and the two energy bands (20 - 50 keV of 
SuperAGILE and 2-12 keV of RXTE/ASM), estimating 
the rate with the time resolution of one orbit 3 ks of 
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net exposure) for SuperAGILE and one dwell 90 s) for 
RXTE/ASM. The SFs of the SuperAGILE data are shown 
in Fig. m from top (OBI) to bottom (0B4). Similarly, the 
structure function of RXTE/ASM data is shown in Fig.fTUl 
with the same disposition of OBI on top and 0B4 at the 
bottom. 

The first evident feature in the structure function 
in both energy bands is the absence of the "horizontal 
branch", indicating that the variability of Cyg X-1 is cor- 
related even at lags of the order of 10^ s, the smallest 
timescale sampled in our observations. The time lag val- 
ues in the data span more than two orders of magnitude, 
from ~ 10^ s up to ^ 4 X 10® s, with the 0B2 the only ex- 
ception, limited to ~ 2 x 10® s, since this observing period 
is shorter than the others. 

The structure functions in the hard X-ray band show 
a monotonic increase between time lags of about 10^ up 
to ~ 10® s. Above ~ 10® s the curve is less smooth, many 
fluctuations appear and, in some cases, a few minima are 
found. Similar behaviour is qualitatively found in the soft 
X-ray band, and also in this case the monotonic region ex- 
tends between 10^ s and ~ 2 x 10® s. In the OBI and 
0B2 the structure function in the two energy bands are 
comparable, while in the 0B3 and 0B4 the SuperAGILE 
data are different from the RXTE/ASM data, which show 
a flatter power-law, with less pronounced fluctuations. 

We define the minima as these values for which the 
structure function is more than 3cr far from the smoothed 
value. With this method we find minima only in the 0B2, 
at 1.42 x 10® s (corresponding to 16.5 days) in SuperAGILE 
and at 0.57 x 10® s (6.7 days), 1.15 x 10® s (13.3 days), and 
1.53 X 10® s (17.7 days) in RXTE/ASM. Interestingly, these 
three minima are all multiples of the same At of 6.7 days, 
of which 13.3 days is the double and 17.7 days has a ratio of 
4/3. The minimum in the SuperAGILE data has a ratio of 
5/2 with the At of 6.7 days. We do not find minima in the 
other observing periods, only fluctuations near the time cor- 
responding to half of the duration of the observation, where 
the scattering of the data increases. This is the case, for ex- 
ample, of the broad "valleys" in the SuperAGILE SFs of 
the 0B3 (around gap of 4 x 10® s) and 0B4 (about 3 x 10® 
s) or the wide fluctuation in the RXTE/ASM sample of the 
OBI (between 2 and 3 xlO® s). 

To confirm our results in hard X-rays, we repeated 
the analysis of the structure function using the publicly 
available Swift/BAT data, in an energy band similar to 
SuperAGILE (15 - 50 keV). Searching with the same sig- 
nificance level, we found minima at 0.42 x 10® s (4.9 days), 
0.99 X 10® s (11.5 days), and 1.42 x 10® s (16.5 days) in 
the 0B2. The last (1.42 x 10® s) is at the same position 
as the one found in the SuperAGILE data. The analysis of 
the Swift/BAT data corresponding to the OBI, 0B3, and 
0B4 confirms the absence of minima in the SF. 

To verify that the structure function is indeed insen- 
sitive to windowing and aliasing of the time series, (e. g. 
^ughes et al. 1992), we computed the SF of a simulated 
time series. At the same time values of the 0B2, shown 
in Fig. [3j we simulated the flux by extracting with a ran- 
dom uniform distribution between the minimum and max- 
imum measured values. Then, we computed the structure 
function of this simulated time series and applied the same 
algorithm used to find the position of the minima. After 
10"* simulations, our algorithm found a fake minimum only 




Fig. 9. Superposition of the first-order structure function 
of the SuperAGILE flux extracted from the orbital integra- 
tion. From top OBI to 0B4 are shown. The arrow mark 
the position of the minimum 



in 71 occurrences, which is the 0.7 %, thus confirming the 
goodness of our results. 

Finally, we fitted the function log[£'^(T)], independently 
in the four observing blocks and the two energy bands, to 
estimate the Hurst exponent H. Only in two cases, 0B3 in 
the SuperAGILE data and OBI in the RXTE/ASM data, 
the structure function shows a break and it is not possible 
to apply the same fit to the complete curve, thus the in- 
terval has been split in two, fitted separately. In the other 
cases, the structure function is described by a single fit. 
The ranges of time lag in which the fits are applied and the 
resulting value of the H exponent are reported in Table |3] 
(SuperAGILE data) and Tablel (RXTE/ASM data). In all 
the observing blocks and in both energy bands, H is always 
significantly lower than the critical value of 0.5, indicating 
antipersistence in the source emission. 

4.4. Alternative methods of timing analysis 

We adopted additional and alternative methods of timing 
analysis to verify the timescales of the system found in the 
structure function: Lomb-Scargle periodogram and auto- 
correlation function. These tools were chosen to comply 
with the structure of our time series, which is very unevenly 
spaced. 

The Lomb-Scargle periodogram is an algorithm specif- 
ically developed to calculate the power sp e ctrum of such 
non-uniformly spaced d ata (jLombl Il976t IScargld 119821 : 
iHorne fc Baliunai |1986[ ). This tool has been already 
adopted in the analysis of the long -ter m variability of 
CygX-1, for example by iRicd (|2008h and iBenlloch et"aI1 
(200^ , to search for the superorbital period. Only some 
of the timescales from the structure function are also 
present in the periodogram of the 0B2: 6.7 and 13.3 days 
in the RXTE/ASM data and 4.9 and 16.5 days in the 
Swift/BAT data. Other timescales, in contrast, are not 
found in the Lomb-Scargle diagram, such as the 16.5 days 
in SuperAGILE, the 17.7 days in RXTE/ASM and the 11.5 
days in Swift/BAT. 
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Fig. 11. Sum of the GRID count maps of the four observ- 
ing periods, The data are smoothed with a Gaussian shape 
of 3 pixel standard deviation, and the circle at the posi- 
tion of Cyg X-1 has a radius of 1.5°, corresponding to the 
instrument point spread function at ^ 200 MeV. 



Fig. 10. Superposition of the first-order structure function 
of the RXTE/ASM flux extracted from the "dwells" of 90 
s duration. From top OBI to 0B4 are shown. The uncer- 
tainty on the SF is generally smaller than the symbol used 
in the plot. The arrows mark the position of the minima. 



By using the Lomb-Scargle periodogram, iRicd (120081^ 
finds a superorbital modulation in the flux of Cyg X-1 with 
a period of 326 ± 2 days. The SuperAGILE data span a 
total duration of ~ 395 days, but the observing periods are 
shorter and not evenly distributed during this time. For 
this reason we cannot significantly test the presence of such 
superorbital modulation, as confirmed by our attempts to 
fold our data on the periods mentioned above. 

We calculated the autocorrelation function from the 
SuperAGILE and RXTE/ASM data in the four observing 
blocks (defined in Table[T]) after removing the dips produced 
by the inferior conjunction of the binary syste m and after 
rebinning the data at six hours. Similar to .Maccarone et al.l 
(|2000D , we computed the autocorrelation function in three 
energy bands: 20 - 50 keV of SuperAGILE, 1.5 - 5 keV 
given by the sum of the A and B bands of RXTE/ASM 
and 5-12 keV from the ASM C band. Apart from the 
peak at zero lag, no other peak is found, indicating that 
the source emission does not show periodicity in all energy 
ranges. Some differences in the width of the autocorrelation 
are found depending on the energy band, but our analysis 
does not show any definite trend. 



5. Results in the gamma ray energy band 

We analysed the GRID data to search for emission of Cyg 
X-1 above 100 MeV using different timescales: one day, five 
days, a single observing block, and the complete AGILE 
observation. For the timescales of one and five days, we 
analysed all the contiguous time intervals included in the 
four observing blocks independently. Compared to the real 
time, the live time is typically ^ 40 ks and ~ 200 ks for 
the two types of integrations, respectively, owing to the 
Earth occultation and the satellite passages through the 
SAA. Cyg X-1 is not detected as a point-like source on any 
timescale and the 2a upper limit depends on the duration 
of the observation: ~ 100 x 10~® ph cm~^ in one day, 
and 50 X 10"® ph cm"^ s"^ and ~ 30 x 10~® ph cm"^ s"^ 
in one observing block. The value of the upper limit also 



depends on the position of the source within the GRID 
field of view (see Table [1]) , which influences the effective 
area, particularly in the last part of the OBI, where the 
upper limits are ~ 30 % higher. 

To find our deepest upper limit on the source flux, we 
accumulated the maps of the whole observation, as shown 
in Fig.[TT] The circle at the position of Cyg X-1 has a radius 
of 1.5°, corresponding to the GRID point spread function 
at ~ 200 MeV. The source is not detected and the upper 
limit above 100 MeV is ~ 5 x 10~* ph cm~^ s~^. From the 
gamma ray map, it is possible to see that the closest source 
is PSR J2021-h3651, located at a distance of 4.6° from Cyg 
X-1, well outside the point spread function. We can then 
exclude the contamination by field sources. 



6. Discussion 

The AGILE mission observed Cyg X-1 4 month- long 
continuous periods of time during the science verification 
phase and the observing cycle 1, between July 2007 and 
December 2008, for a total net exposure to the source of 
about 6 Ms. Because of the AGILE pointed strategy, dur- 
ing these observing periods the monitoring of the source is 
continuous, interrupted only by the periodic occultation by 
the Earth. Under this condition, the SuperAGILE data on 
Cyg X-1 are a unique set in the hard X-rays. Other hard 
X-ray experiments provided monitoring data on Cyg X-1, 
but all of them are either (i) sparse, short snap s hots of the 
sourc e e.g. C GRO/BATSE ("Crarv et al.l [1991 iLing et al.l 
119971 : iBrock sopp et al JJ99; Zdziars kreTalT 120021 ) and 
Swift/BAT (Rico 200j) or (ii) short pointed observations 
(e.g., INTEGRAL/ISGRI: [Cadolle Bel et al.l [20061) . By 
taking as a reference a colour-colour diagram built with the 
data of the RXTE/ASM in the soft X-ray energy range, we 
verified that Cyg X-1 remained in its low-hard state the 
whole time. 



6.1. X-ray variability on short timescales 

We used the SuperAGILE data to systematically inves- 
tigate the source variability over our 160 days of mon- 
itoring on different timescales, from ~ 300 s to the ~ 
1-hour timescale, up to week-long time segments. For 
the time periods when the source exhibited a signifi- 
cant variability on the minute timescale, we also (un- 
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successfully) searched for potential flares on the tens-of- 
second to subsecond-timescale in the SuperAGILE raw 
lig htcurves. Such variabili t y has been observed in the past 
by iGierlinski fc Zdziarskil (|2003D on the ms timescale in 
the soft X-rays with RXTE/PCA, and by Golenetskii et 
al. (2003) on the hour timescale in the hard-X/gamma 
rays using the experiments of the Interplanetary Network. 
Both types of variability were observed with the source 
in both soft and hard spectral states in a flux range of 
a few X 10"'' erg cm~^ (corresponding to a luminos- 
ity of ^ 10'^* ergs""'^, assuming a distance of 2 kpc) in 
the energy ranges 3-30 keV and 15 - 300 keV, respec- 
tively. During our observations we detected several flares 
on timescales of a few hours to half a day, at times corre- 
sponding to MJD 54420.4, 54585.6, 54604.6, 54607.9 and 
54646.5. Contrary to the cases cited above, which involved 
flux variability of one order of magnitude, these flares cor- 
respond to flux increases of approximately a factor of 2 - 
3, passing from ~ 5 — 6 x 10^^ erg cm^^ s~^ to maxima of 
about ^ 1.4 X 10~* erg cm~^ s~^, as measured on the hour 
timescale in the 20 - 50 keV energy range. Results si milar 
to ours were also reported bv iBrocksopp et aD (|1999[ ) in a 
study based on dat a of BATSE and RXTE/ASM, and by 
iMalzac et ahl (pOOSi ) using INTEGRAL/ISGRI (in the en- 
ergy band 18 - 40 keV). In particular, th e 18 May 2008 
MJD 54604) flaring episode reported by iNeronov et al.l 
20081) . with a duration of about 8 hours, is among the 
flares observed by SuperAGILE during the 0B3. The flux 
measured by SuperAGILE in 20 - 50 keV between 15:02 
and 15:36 UT, when INTEGRAL/ISGRI detected the max- 
imum of the source emission, is ~ 1.5 Crab. 

6.2. Long-term variability 

Interestingly, with the exception of the flares discussed 
in the previous section, during the SuperAGILE moni- 
toring the hard X-ray flux of Cyg X-1 remained quite 
stable around ^ 5 — 6 x 10~^ erg cm~^ s~^, with fluc- 
tuations on much longer timescales. We investigated this 
long-term variability by using different analysis tools. The 
SuperAGILE data are characterized by long gaps between 
the observing periods, and shorter gaps due to the Earth's 
occultation in every orbit. For this reason we applied the 
analysis of the structure function for the first time in the 
study of a Galactic binary source, but frequently used for 
the active galactic nuclei in th e radio and gamma ray energy 
bands (see for exampl e Hughe s~et al.lll992l : ISimonetti et al.l 
119851: iNandikot kur et al.l il9971. We performed a compara- 
tive analysis in the hard and soft X-rays, complementing 
the SuperAGILE data with those of RXTE/ASM available 
in the public archive. Since the slope of the structure func- 
tion is connected to the probability of long duration trends 
in the emission mechanism, while the minima in the same 
function are related to the typical timescales of the same 
process, we address the two topics separately. 

By using a logarithmic fit to the structure function ver- 
sus the time lag, we find that, throughout the whole period 
of our analysis, the Hurst exponent H is significantly be- 
low the critical value of 0.5, indicating that the source is 
antipersistent. This antipersistence is a property of random 
non-stationary processes and indicates that an increase in 
the p ast is likely t o be followed by a decrease in the fu- 
ture (jKarneij |2005[ ). Antipersistent systems are not com- 
mon in nature and show a dominant negative feedback at 



work that produces "oscillations" , i. e. repetition of similar 
features but without periodicity, because the distance in 
time is no t const ant ( see for exampl e Koutsoyiannis 200^. 
iGil-Alanal (|2006D and lKarneil (|2005[ ). for example, find an- 
tipersistence in the climatological study of the time series 
of the temperature in various layers of the Earth atmo- 
sphere and discuss possible mechanisms of negative feed- 
back ( ice albedo, water vapour and clouds). On the other 
hand, lAlvarez-Ramirez et al.l (|2002[ ). always by means of 
the Hurst exponent, find that the crude oil market is a per- 
sistent process with long-run memory effects and different 
timescales at work: the "oscillations" from days to weeks 
are apparently generated by the action of market specu- 
lators and are superimposed on a long-time trend, with 
characteristic timescale of weeks to three months. 

The description of the antipersistent systems matches 
the behaviour of Cyg X-1, whose luminosity does not seem 
to vary on long timescales and shows non-periodic "oscilla- 
tions" around an average value of ^ 2 x 10^^ erg s~^. This 
aperiodic variation is well known especi ally in the low / hard 
spectral states, (see for example Makish ima et al.|l2008l and 
references therein). A complete discussion of the possible 
negative feedback mechanisms in the emission of Cyg X-1 
would require a detailed analysis of the interactions be- 
tween the various components at work (stellar wind, accre- 
tion disk, Comptonizing plasma, and relativistic jet) and is 
beyond the scope of the present paper. 

The position of the minima in the structure function 
i s an indication of the typical timescales of the source 
(INandikotkur et aI1ll997f) . In our analysis we only found 
minima in the structure function of the observing block 
2 (Aprfl 2008), but at different positions in SuperAGILE 
(16.5 days) and RXTE/ASM (6.7, 13.3, and 17.7 days). 
The analysis of the structure function of the Swift/BAT 
data in the same periods confirms the SuperAGILE mini- 
mum at 16.5 days and is accompanied by two more min- 
ima on slightly different timescales (4.9 and 11.5 days). 
The minima that we detected do not correspond to any 
known timesca le of Cyg X-1, whi ch has an orbital period 
of ~ 5.6 days (jLasala et al.l 1199 8^. and are not confirmed 
by other techniques of timing analysis. In fact, we do not 
detect any periodicity using the autocorrelation function, 
and we found only some of the timescales in the Lomb- 
Scargle periodogram. An inspection of the lightcurves of 
the 0B2 reveals that the timescales found in the structure 
function correspond to features in the time series: 6.7 days, 
which is also present in the Lomb-Scargle periodogram, is 
the average distance between the minima in the soft X- 
ray lightcurve and 16.5 days is the separation of the two 
maxima of the emission in hard X-rays. From all these re- 
sults we conclude that the timescales that we found in the 
structure function correspond to emission features of the 
specific observation period and not to general properties of 
the source. 

6.3. The searcli for gamma ray emission 

Thanks to the long and continuous monitoring of AGILE, 
we were able to investigate the possible emission of gamma 
rays from Cyg X-1 on both short and long timescales. An 
investigation on the daily timescale did not provide any 
significant detection. We then increased the exposure time 
to five days, one observing block (~ 1 month) and finally 
the whole 160-days observation. We did not find any signifi- 
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cant emission of gamma rays on any of these timescales, the 
typical 2(7 upper limit being ~ 100 x 10^^ ph cm~^ s~^ on 
the one-day and ~ 50 x 10~^ ph cm~^ s^^ on the five-day 
timescale. The tightest 2a upper limit derives from the in- 
tegration over the whole observation (6T Ms net exposure) 
and corresponds to ^ 5 x 10^^ ph cm^^ s^^. 

Recently, an episodic flare from Cyg X-1 was detected 
in gamma rays by AGILE in an interval of time subse- 
que nt to that consi d ered i n the present paper, and reported 
by ISabatini et al.l (j201Clt ). Between 15 and 16 October 
2009 (MJD 55119.97 - 55120.96), a flux of (232 ± 66) x 
10~^ ph cm~^ s~^ above 100 MeV was measured, at a sig- 
nificance level of 5.3a pre-trial and 4 a post-trial, when the 
source was in a hard spectral state. ISabatini et al.l (|2010[ ) 
also extended the search for persistent emission over all the 
AGILE archival data of the Cygnus Field, including the in- 
terval covered by the present paper, and did not find any 
significant detection. Given the longer integration time, the 
authors also obtained an upper limit slightly deeper than 
found in our analysis, ~ 3 x 10~^ ph cm~^ s~^, with an 
integration time of about 300 days, spread over 2.5 years. 

Cyg X-1 has also been detected at TeV energies by 
MAGIC ([Albert et al.ll2007D during an ~hour flare, at the 
rising edge of a flaring episode detected in both soft and 
hard X-rays. This rare, and so far unique, event occurred 
with the source in a hard spectral state, as in our observa- 
tions, although reaching a flux of '^2 Crab in hard X-rays, 
while our X-ray data never show the source above 1.5 Crab. 
The MAGIC observation shows a power-law spectrum be- 
tween 150 GeV and 2 TeV with photon index 3.2 ± 0.6, 
yielding a flux of 7 x 10~^^ ph cm~^ . If we assume 
that the spectrum can be extrapolated to the energy range 
0.1 - 1 GeV, as when the emission site is far from the 
disk and the plasma, so that the cross-section of the pho- 
ton conversion into e lectron positron pairs is negligible (see 
IZdziarski et"alll2009l ). we obtain an expected flux value of 
~ 63000 x 10"^ ph cm-2 8"^ Even "diluting" such a flux 
on a 24-hour timescale, it would correspond to an expected 
day-averaged flux of ^ 3000 x 10"^ ph cm^'^ s^^ , well above 
the GRID daily sensitivity of ^ 100 x 10"^ ph cm"^ s~^ 
Thus, regardless of the spectral state of the source, our anal- 
ysis can exclude the presence of gamma ray (GeV-)flares 
similar to what is observed by MAGIC at TeV energies, 
during the periods of the AGILE observation. 

Cyg X-1 was observed by the instruments aboard the 
Compton Gamma Ray Observatory on several occasions 
during the flrst three cycles (1 991 - 1994), when th e 
source was m a low/hard state (jMcConnell et all l2000l) . 
COMPTEL measured an energy spectrum described by a 
simple power-law with photon index 3.3 ± 0.4, between 
0.75 and 5 MeV. Above 5 MeV, the statistical quality 
of the data does not allow placing more than an upper 
limit consistent with the extrapolation of the spectrum 
at lower energies. During the same observations, EGRET 
did not detect the source, with an upper limit on the flux 
at 10 X 10-8 ph cm -^ s"^ (100 - 200 M eV). The deepest 
AGILE upper limit (jSabatini et al.|[2010l) is about a factor 
of three lower than EGRET. However, an extrapolation of 
the COMPTEL spectrum above 100 MeV would provide an 
expected flux of 5 x 10^^ ph cm^^ s^^, quite consistent 
with both of them. Our data, then, do not provide evidence 
of any spectral cut-off above the COMPTEL energy range. 

Interestingly, Cyg X-1 was also detected by COMPTEL 
in high/soft state, up to about 10 MeV energy 



(jMcConnell et al.ll2002D . At that time the EGRET instru- 
ment was switched off and no GeV observation is available. 
The COMPTEL spectrum of the soft state can be mod- 
elled with a single power-law of photon index 2.6 extend- 
ing from 1 to 10 MeV. Extrapolating this energy spectrum 
to the AGILE energy band would provide a flux expecta- 
tion of 30 x 10"^ ph cm~^ s~^, easily detectable by the 
AGILE/GRID on a timescale of a few days. This can be 
seen as an additional indication that Cyg X-1 most likely 
remained in its low/hard spectral state during all the time 
of the AGILE. 

The AGILE upper limits on GeV-emission from Cyg X- 
1 range from 100x10^* ph cm~^ (one day timescale) 
to 3x10^8 ph cm^^ s^^ (^ 2.5 years average ). Assuming 
a source distance of 2 kpc, (jZioIkowskil |2005| ) they corre- 
spond to an isotropic source luminosity of ~ 8 x 10"^^ erg 
s~^ and '-^ 2 X lO^'^ erg s""'^, respectively. Different authors 
elab orated models for the high-energy emission from Cyg 
X-1. IZdziarski et al.l (j2009t ) proposed a model t o explain the 
short TeV flare detected in 2006 by MAGIC ([Albert et al.l 
2007). They propose that high-energy photons are gener- 
ated by electrons accelerated to TeV energies close to the 
central X-ray source. The interaction of such TeV photons 
with the photon fleld of the X-ray source yields produces 
e^ pairs. They find that, for initial photon energies above 
~3 TeV, photons can travel far enough, initiating extended 
pair cascades, in turn producing an observable flux. The 
predicted photon energy spectrum depends on the energy 
of the primary electron injection and on the stellar tem- 
perature. Under the conditions observed in X-rays during 
the TeV flare (when the hard X-ray flux was only ~^two 
times higher than what we observed with SuperAGILE, 
iMalzac et aI1l2008D . the model predicts a flux in the AGILE 
energy range as high as 4 x 10'^'* erg s^^. This is a lumi- 
nosity value higher than the AGILE time-averaged upper 
limit on the persistent emission, while it is still compatible 
(2 times less) with the 1-day upper limit. 

Alternativ e mod els have also been proposed. 
lAraudo et al.l (|2009f) studied the interaction of the 
jet of Cyg X-1 with the wind of the companion, expected 
to be clumpy, and predicted the possible spectral energy 
distributions under different conditions of the system. Since 
the radiation is considered to be produced by particles 
accelerated in the shock due to the interaction between the 
jet and a cloud, the emission is predicted to be transient, 
with the duration timescale expected in the range between 
minutes to a few hours. The frequency of occurrence of 
such transient emission depends on the clumps' number 
and size, up to appearing as a modulated steady emission. 
Interestingly, under some magnetic field / clump size 
conditions, the model predicts the radiative output due 
to inverse Compton effect to reach luminosities in the 
AGILE energy range as high as ~ few x lO'^* erg s~^, 
comparable to the AGILE daily upper limits and nearly 
one order of magnitude h i gher than the yearly-average 
limit. Also lOreUana et al.l |2007t) suggest that the GeV 
- TeV emission in microquasars may originate from the 
interaction between the jet and the stellar wind of the 
companion. In this case the inelastic interactions between 
the relativistic protons in the jet and the cold protons 
of the stellar wind produce charged (tt^) and neutral 
(tt'') pious. The jet is assumed to be continuous and to 
contain a randomly oriented magnetic field. Two main 
channels of gamma ray production emerge from the model: 
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the decay of neutral pions and the emission by charged 
leptons, resulting from the decay of charged pions and 
from the photon-photon interactions, producing photons 
via synchrotron emission and inverse Comp ton scattering 
on the l ow energy stell a r pho tons. Similar to lAraudo et al.l 
(|2009() . lOrellana et al.l ()2007() present the spectral energy 
distribution (SED) expected in their model for the two 
emission channels (tt^ and 7r°) using the parameters of Cyg 
X-1. The gamma ray luminosity produced by the decay of 
the neutral pions depends on the subtended solid angle 
and on the inclination of the jet toward the line of sight. 
Taking an inclination of 30° into account (measured by 
iGallo et "aI1l2005,. in radio) and a semi-opening angle of 0.1 
radians, the expected luminosity is more than three orders 
of magnitude lower than the AGILE upper limit, at an 
average energy of 100 MeV. On the other hand, the leptons 
emitted in the decay of the charged pions can produce 
gamma rays via inverse Compton scattering of low energy 
photons from the companion star. The interactions of 
gamma rays with low energy photons may in turn produce 
relativistic pairs that can Compton upscatter the low 
energy photons or trigger an electromagnetic cascade. In 
the SED of such secondary emission from charged pions, 
adopting two values of the magnetic field, 10^ G and 
10^ G, the luminosity at 100 MeV is slightly above the 
^ 2 X 10^-^ erg s~^ AGILE upper limit for the highest 
magnetic field (10^ G) while the model is still compatible 
with our data for a lower magnetic field, 10"* G, whose 
luminosity is about one order of magnitude less than our 
upper limit. 

7. Summary and conclusions 

We reported the first campaign of observation of the hard 
spectral state of Cyg X-1 in gamma rays (> 100 MeV) on 
various timescales, with exposures ranging from one day 
up to ~ 160 days. We monitored the source simultane- 
ously in hard (20 - 50 keV) and soft (2 - 12) X-rays using 
SuperAGILE and RXTE/ASM. The observation in hard 
X-rays by SuperAGILE shows the well-known erratic vari- 
ability of the flux of Cyg X-1. The analysis of both the 
hardness ratio, estimated from the SuperAGILE flux in the 
two energy bands 20 - 25 keV and 25 - 50 keV, and of the 
colour-colour diagram, obtained from the pub lic data of 
RXTE /ASM following the method reported by iReig et al.l 
(|2002tl . does not show any transition of spectral state. We 
adopted exposure times ranging from ^ 300 s up to one 
day and found the typical short time flickering of Cyg X-1 
with intensity variations by about a factor of two. 

From the study of the first-order structure function, we 
did not find any short lag plateau ("horizontal branch"). 
A power-law behaviour is found, with the Hurst exponent 
smaller than 0.5, denoting that the emission mechanism 
of Cyg X-1 is antipersistent, i. e. dominated by negative 
feedback. We also found timescales from the minima in the 
structure function but, from the analysis of the autocor- 
relation function and the Lomb-Scargle periodogram, we 
derive that these timescales are more probably related to 
particular patterns in the specific lightcurve, such as the 
distance between the repeated minima or the peaks of the 
emission, rather than to intrinsic properties of the source. 

Cyg X-1 is not detected as a point like source above 
100 MeV, and we find values of the upper limit rang- 
ing from ^ 100 X 10~^ ph cm~^ s~^ in one day down to 



^ 5 X lO^'* ph cm^^ s^^ for the whole observation, about 
a factor of two lower than from the EGRET data. We 
compared the luminosity derived from the AGILE upper 
limit above 100 MeV (assuming a distance of 2 kpc to 
the source) with various models of GeV - TeV emission. 
The predictions of th e pairs' cascade model (proposed by 
'Zdziars ki et al.l [20091 to explain the TeV flaring emission 
detected by MAGIC) are compatible with our upper limit 
only for flaring emission (e. g. one-day timescale), while 
they are about one order of magnitude higher when com- 
pared to the year-long upper limit. We compared our results 
with the predictions of two alternative models, involving 
hadronic interactions between the cold matt er of the stel- 
lar wi nd an d the relativistic jet, p roposed bv lAraudo et al.l 
()2009[ ) and lOrellana et al.l (|2007l) . The first model is more 
suitable for flaring emission, while the second one can ex- 
plain the persistent emission better. We find that the lu- 
minosity in the ~ 100 MeV range is compatible with the 
results of the AGILE monitoring for either short transient 
or for selected model parameters. 
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Obs. 


Date 


Modified 


off-axis angle 


exposure 


period 






Julian Day 


[degrees] 


[ks] 


OBI 


02 Nov 2007 13:50 UT - 


16 Dec 2007 10:27 UT 


54406.576 - 54450.435 


2.4 - 27.6 


1634 


0B2 


10 Apr 2008 15:10 UT - 


30 Apr 2008 11:15 UT 


54566.632 - 54586.469 


18.3 - 15.4 


893 


OB3 


10 May 2008 17:52 UT - 


30 Jun 2008 11:10 UT 


54596.744 - 54647.465 


4.6 - 32.7 


1934 


OB4 


17 Oct 2008 17:00 UT - 


01 Dec 2008 11:37 UT 


54756.708 - 54801.484 


15.1 - 18.0 


1632 





Fig. 3. SuperAGILE lightcurve of Cyg X-1 during the observing periods 1 (top left), 2 (top right), 3 (bottom left) and 

4 (bottom right). Each point represents the normalized rate estimated in one satellite orbit (with exposure of 3 ks). 
The normalized rate values arc converted into units of erg cm^^ s^^ considering the source average spectrum. Only the 
statistical uncertainty is shown on the points in the plots. The average systematic uncertainty on the flux calibration is 
plotted in the top right angle of each plot. 

Table 2. Parameters of the spectral fit of Cyg X-1. Xr is the reduced chi square, the uncertainties are at 90 % confidence 
level and include a 5 % systcjniatic cTror. 



Normalized rate 


photon index 


flux 


xi 


degrees of freedom 


[cts cm^'^ s^"'] 




[erg cm^^ s^^] 






0.12 


1 63+""'* 


(7.6 ±0.3) X 10"" 


1.159 


31 


0.17 


1.60 ± 0.11 


(8.3 ±0.4) X 10"^ 


1.113 


32 


0.08 


1.74 ± 0.09 


(4.0 ± 0.4) X 10 


0.986 


29 
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Fig. 6. SuperAGILE daily lightcurve of Cyg X-1 during the observing blocks 1 (top left), 2 (top right), 3 (bottom left), 
and 4 (bottom right). Each point represents the normalized rate estimated in one day (with exposure of ~ 40 ks). The 
rate values arc converted into units of erg cm~^ s~^ considering the source average spectrum. The hardness ratio is 
computed as the ratio of the emission in the HE band (25 - 50 keV) to the LE one (20 - 25 keV). Only the statistical 
uncertainty is shown on the points in the plots. The average systematic uncertainty is plotted in the top right corner of 
each plot. 



Table 3. Analysis of the structure function in the hard X-ray band 



Obs. start lag end lag H 

period [s] [s] 



OBI 


1Q4.9 


10" 


0.185 


OB2 


1Q3.9 


10.5.2 


0.147 


OBS 


1Q3.9 


104.9 


0.115 




1Q6.2 


IOS.6 


0.210 


OB4 


104.9 


105.7 


0.170 



Table 4. Analysis of the structure function in the soft X-ray band 



Obs. 


start lag 


end lag 


H 


period 




[s] 




OBI 


103.V 


10^'^ 


0.144 




105.2 


10" 


0.186 


0B2 


10-* 


106.2 


0.148 


OBS 


lO^-^ 


105.2 


0.066 


OB4 


10*-2 


105.2 


0.050 
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Fig. 7. Superposition of the Super AGILE (top panels) and RXTE/ASM (middle panels) lightcurves of Cyg X-1 and 
hardness ratio (bottom panels) during the observing periods 1 (top left), 2 (top right), 3 (bottom left), and 4 (bottom 
right). The arrows show the position of periods of increase in the source emission that are visible in the hard X-ray band 
but not in the soft one. Both lightcurves in each period are accumulated on the same time interval with a bin size of six 
hours. Only the statistical uncertainty is shown on the points in the plots of the SuperAGILE flux and of the hardness 
ratio. The average systematic uncertainty is plotted in the top right corner of each plot. 



